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Abstract 

A two-body system hypothetically affected by an additional radial acceleration Hv r , 
where v r is the radial velocity of the binary's proper orbital motion, would experience long- 
term temporal changes of both its semimajor axis a and the eccentricity e qualitatively 
different from any other standard competing effect for them. Contrary to what one might 
reasonably expect, the analytical expressions of such rates do not vanish in the limit M — > 0, 
where M is the mass of the primary, being independent of it. This is a general requirement 
that any potentially viable physical mechanism able to provide such a putative accelera- 
tion should meet. Nonetheless, if H had the same value Ho of the Hubble parameter at 
present epoch, such rates of change would have magnitude close to the present-day level 
of accuracy in determining planetary orbital motions in our Solar System. A tension with 
recent observations may even be present for Mercury and Mars. However, general relativ- 
ity, applied to a localized gravitationally bound binary system immersed in an expanding 
Fricdmann-Lemaitre-Robertson- Walker, does not predict the existence of such a putative 
radial acceleration at Newtonian level. Instead, it was recently shown in literature that 
an acceleration of order H and directed along the velocity v of the test particle occurs at 
post-Newtonian level. We worked out its orbital effects finding well-behaved secular rates of 
change for both a and e proportional to the Schwarzschild radius r s of the primary. Their 
magnitude is quite small: the rate of change of a amounts to just 20 [ira per century in our 
Solar System. Finally, we discussed certain basic criteria of viability that modified models 
of gravity should generally meet when their observable effects are calculated. 

PACS: 04.80.-y; 04.80.Cc; 98.80.Es; 95.10.Km 

1 Introduction 

In this paper, we first deal with a certain hypothetical anomalous radial acceleration propor- 
tional to the radial velocity of the orbital motion of a two-body system through a coefficient H 
having dimensions of T -1 . 

Some intuitive, Newtonian-like guesses about the possibility that such a putative acceler- 
ation may exist as a local manifestation of the cosmic expansion in the case of non-circular 
motions are offered in Section [2j They are motivated by the known fact that, within certain 
limits, several well established key features of a homogeneous and isotropic expanding universe, 
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and also of its influence on local gravitationally bound systems, can be practically inferred 
within a classical framework [H [21 [3j [H [5] . 

Such a putative extra-acceleration may, in principle, have interesting phenomenological 
consequences since it would induce peculiar orbital signatures which would not be mimicked by 
any other known competing dynamical effect. Moreover, by assuming for H a value equal to 
that of the Hubble parameter at present epoch, the magnitude of these exotic effects for the 
planets of our Solar System would be close to the current level of accuracy in determining their 
orbits. Tensions might even occur between data and predictions in the case of Mercury and 
Mars. These topics are treated in Section [3] 

In Section H] we try-unsuccessfully-to find a theoretical justification for the guesses of Sec- 
tion [2] rooted in a full general relativistic treatment of the orbital dynamics of a local system 
embedded in an expanding homogeneous and isotropic Friedmann-Lemaitre-Robertson- Walker 
(FLRW) spacetime metric. Nonetheless, our results remain valid from a phenomenological 
point of view because of their actual independence of any specific theoretical scheme, and can 
be viewed as observational constraints on such a putative exotic force, whatever the physical 
mechanism yielding it (if any) may be. Moreover, we feel that the numerical values coming out 
from our analysis are interesting if compared with the observations, and may pursue further 
investigations to find a possible physical origin, cosmological or not, for it. Cosmological effects 
linear in H were recently derived by Kopeikin [6] for the propagation of electromagnetic waves 
between atomic clocks in geodesic motion in a FLRW background. 

In Section [5] we build on the certain aspects discussed in previous Sections and provide some 
very general viability criteria that must be met by modified models of gravity [7] in order not to 
give rise to unphysical observable effects. In particular, it is important to check the behaviour of 
their detectable predictions in the limits G — > 0, M — > 0, where G is the Newtonian constant of 
gravitation and M is the mass of the central body acting as localized source of the gravitational 
field. 

Section [6] provides an overview of the results obtained. 

2 A local cosmological effect linear in H for two-body orbital 
dynamics? 

Let us define the Hubble parameter H and the deceleration parameter q in the usual way as 



where S(t) is the cosmological scale factor. From an observational point of view, the value of 
the Hubble parameter at present epoch is [8] 




(1) 




(2) 



H = 100 km s" 1 Mpc _1 /i = (73.8 ± 2.4) km s _1 Mpc 



-i 



(3) 



so that 



h = 0.738 ± 0.024. 



(4) 
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The deceleration parameter can be connected to directly determined quantities from observa- 
tions by means of [9] 

qo = —5 ^A,o, (5) 

where f2 m o and Oa,o are the current matter and dark energy densities, respectively, in units of 
the critical density [9] 

crit _ 3 H ( ^ 

From [TO] 

n mfi h* = 0.13341°;™!, (7) 



n A) o = 0.728±g;gii, (8) 

and eq. Q, it turns out 

q ~ -0.7. (9) 

Let us now consider a two-body system composed by a central object of mass M and a 
test particle gravitationally bound to it. Their proper motion is determined by their mutual 
gravitational interaction according to the Newtonian/Einsteinian laws. It superimposes on the 
global Hubble flow in such a way that the radial velocity of the test particle is the sum of such 
two contributions 

v r =r = ^ orb ) + Hr. (10) 

In it, r is their relative distance and 

v {orb) = n b ae sin / 



e 



is the radial component of the velocity vector v of the test particle along its standard two-body 
Keplerian ellipse where a is the semimajor axis, e is the eccentricity, / is the true anomaly, and 
nb = \J GM/aP is the Keplerian mean motion. As shown by Table [H the contribution of the 
Hubble flow to the planetary radial velocities is quite negligible in the Solar System. 

When the accelerations are computed, neglecting the proper motion yields the well known 
Hooke-like term quadratic in H. Indeed, starting from the Hubble law 

f = Hr, (12) 

yields 

A r = Hr + Hr. (13) 

By recalling that 



H = | - H 2 (14) 



and by using eq. §Z§ and eq. (fTZj) one gets just 

^ 2N ) = - qH 2 r. (15) 
If we include vi° rh ^ in the second term of eq. (|13p . we get an additional term linear in H 

= Hv^ (16) 

at Newtonian level. 
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Table 1: Keplerian two-body and cosmological radial velocities of the planets of the Solar 
System, in km s" 1 . As far as the Keplerian two-body velocities are concerned, the maximum 
values of eq. (fTT|) (sin/ = 1) were taken. The average planetary distances d = a (l + ^\ 
were used in the Hubble terms which were evaluated at the present epoch according to [5] 
H Q = 73.8 km s" 1 Mpc" 1 = 2.39 x 10~ 18 s _1 . 



Planet 


Vr (km s l ) 


H d (km s- 1 ) 


Mercury 


10.05 


1 x 10~ 10 


Venus 


0.23 


2 x 10" 10 


Earth 


0.49 


3 x 10" 10 


Mars 


2.26 


5 x 10~ 10 


Jupiter 


0.59 


2 x 10~ 9 


Saturn 


0.52 


3 x 10~ 9 


Uranus 


0.31 


7 x 10~ 9 


Neptune 


0.04 


1 x 10~ 8 


Pluto 


1.22 


1 x 10~ 8 



3 A cosmological test particle acceleration of order Hv r : a phe- 
nomenologically appealing possibility 

Table [2] shows the values of A^ H ^ and Ar ^ compared with the Newtonian monopoles for 
the planets of the Solar System. It can be noticed that the orders of magnitude of Ai HN ^ and 

Table 2: Newtonian and cosmological accelerations of the planets of the Solar System, in m 
s~ 2 . As far as the Hubble terms of order H are concerned, the maximum values of eq. (Ilip 
(sin/ = 1) were taken. The average planetary distances d = a (l + ^\ were used in both 

the Newtonian and the Hubble terms of order H 2 which were evaluated at the present epoch 
according to [5] H = 73.8 km s^ 1 Mpc" 1 = 2.39 x 10" 18 s" 1 and [S] q 0.7. 



Planet 


GM/d 2 (m s~ 2 ; 


1 H vi OTh) (m s- 2 ) 


-qoHfa (m s 2 ) 


Mercury 


4 x 10" 2 


2.40 x 10" 14 


2 x 10" 25 


Venus 


1 x 10~ 2 


5 x 10~ 16 


4 x 10~ 25 


Earth 


6 x 10~ 3 


1.2 x 10~ 15 


5 x 10~ 25 


Mars 


2 x 10" 3 


5.4 x 10" 15 


8 x 10" 25 


Jupiter 


2 x 10" 4 


1.4 x 10~ 15 


2.6 x 10~ 24 


Saturn 


6 x 10~ 5 


1.2 x 10" 15 


4.9 x 10~ 24 


Uranus 


2 x 10" 5 


7 x 10~ 16 


9.8 x 10~ 24 


Neptune 


6 x 10~ 6 


1 x 10~ 16 


1.54 x 10~ 23 


Pluto 


4 x 10~ 6 


2.9 x 10~ 15 


2.09 x 10~ 23 



Ar are completely different, being the terms linear in H of the order of 10 — 10 m 
s~ 2 . Interestingly, such values are neither too large to unrealistically compromise the agreement 
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between theory and observations nor too small to be completely undetectable in any foreseeable 
future. Indeed, according to a recent analysis of various kind of planetary data by Folkner 
the largest unmodelled radial acceleration in the Solar System allowed by observations is just of 
the order of 10 -14 m s -2 . More precisely, from the current Mars radio range data set, Folkner 
found upper bounds on a radial acceleration of Earth and Mars to be less than 3 x 10~ 14 
m s -2 and 8 x 10 -14 m s -2 , respectively. From Cassini radiotechnical data Folkner [11] inferred 
an upper bound of 1 x 10 -14 m s~ 2 for Saturn. Thus, according to Tabled Ai HN ^ seems to be 
not too far from the edge of the present-day detectability. Moreover, the maximum figures for 
A ( r m quoted in Table are not in contrast with the upper bounds by Folkner [11] for Earth, 
Mars and Saturn. On the other hand, it must be remarked that Folkner [11] did not release 
details about the time and/or spatial variability of the anomalous acceleration constrained. 
This fact generally does matter since ^4^™^ is not constant, so that the constraints by Folkner 
[TT] may not be straightforwardly applicable to • 

We now explicitly work out some dynamical orbital effects induced by eq. (116p on some 
quantities routinely determined from observations by astronomers. The standard Gauss equa- 
tions [12] for the variations of the osculating Keplerian orbital elements, applied to eq. (|16p 
evaluated onto the unperturbed Keplerian ellipse by means of eq. (fTTI) . yield the following 
non-zero long-term rates of change of the semimajor axis a and the eccentricity e 

d^ N > = 2aH (l - y/l - e 2 ) = aH e 2 + O (e 4 ) , (17) 



j(flN) 



H {l-e 2 ) (l-VT^ 



H e / 3 



e v ' = : = — — l - -e 



+ O (e 4 ) . (18) 

It is intended that eq. (j!7)l -eq. fjl8j) are averages over one full orbital revolution of the test 
particle. The other orbital elements are left unaffected. Note that, according to eq. (|17j) -eq . 
()18p . both the osculating semimajor axis and the eccentricity increase. This implies that the 
osculating mean distance d = a(l + e 2 /2) 

= 3aH ( 1 - y/T=#) = 3 -^- + O (e 2 ) . (19) 



Rather surprisingly, eq. (|17|) -eq. (|19p . which describe small changes occurring in the otherwise 
close orbit of a gravitationally bound system, do not depend on M: in the limit M —> 0, eq. 
()17p -eq. (|19p do not vanish. The limit G — > does not pose problems in the sense that eq. 
([T7]) -eq. (fl~9j) correctly vanish. Indeed, in a spatially flat FLRW universe, it is [9] 



= *f£, (20) 

where p is the density of the cosmic fluid, inclusive of any dark energy component. In Table 
[3] we calculate eq. (fT7|) -eq. (fT8j) for the planets of the Solar System. It is interesting to 
compare the values in Table [3] with the experimental bounds in Tabled] preliminarily inferred 
from a multi-year fit by Pitjeva [13] for the EPM2006 ephemerides. In general, the predicted 
rates of Table [3] are compatible with the bounds of Table H] although discrepancies occur for 
Mercury and Mars at 4 — a level. It must be stressed that, so far, astronomers did not explicitly 
determine corrections to the standard Newtonian rates of change of a and e from observations: 
the figures in Table U] were inferred rather naively by simply taking the ratios of the formal, 



5 



Table 3: Predicted cosmological rates of change of the semimajor axis a and the eccentricity e 
of the planets of the Solar System for [H] H = (73.8 ± 2.4) km s" 1 Mpc" 1 = (2.39 ± 0.07) x 
10~ 18 s" 1 = (7.5 ± 0.2) x 10~ 9 cty" 1 in eq. (fT7|) -eq. (fTS]) . The rates of a for Mercury and 
Mars, and the rate of e for Mars are larger by one order of magnitude than the corresponding 
experimental uncertainties in Tabled! 



Planet 


• (M) ( m cty -l) 


■{HN) ( cty -i) 


Mercury 


18.68 ±0.61 


(7.5 ±0.2) x 10" 10 


Venus 


0.03 ±0.001 


(2 ±0.08) x 10" 11 


Earth 


0.31 ±0.01 


(6 ±0.2) x 10~ n 


Mars 


15.04 ±0.49 


(3.5 ±0.1) x lO" 10 


Jupiter 


11.94 ±0.38 


(1.7 ±0.05) x 10~ 10 


Saturn 


31.26 ± 1.01 


(2.0 ± 0.06) x 10~ 10 


Uranus 


46.56 ± 1.51 


(1.7 ±0.05) x 10~ 10 


Neptune 


2.47 ±0.08 


(3 ±0.1) x 10" 11 


Pluto 


2849.75 ± 92.67 


(9.0 ±0.3) x 10~ 10 



Table 4: Uncertainties in the rates of change of the semimajor axis a and the eccentricity e of 
the planets of the Solar System. They were inferred by taking the ratios of the formal errors in 
Table 3 of [13], all rescaled by a factor 10, to the data time span AT = 93 yr (1913-2006) of the 
EPM2006 ephemerides used by Pitjeva [13]. The results for Saturn are relatively inaccurate 
with respect to those of the inner planets since radiotechnical data from Cassini were not yet 
processed when Table 3 of [13] was produced. Here cty stands for century. 



Planet 


a a (m cty r ) 


o"e (cty x ) 


Mercury 


3.6 


4 x 10" 9 


Venus 


2.3 


2 x 10~ 10 


Earth 


1.5 


5 x lO" 11 


Mars 


2.8 


5 x 10~ n 


Jupiter 


6612.9 


2 x 10~ 8 


Saturn 


45763.4 


1 x 10~ 7 


Uranus 


433269.0 


3 x 10~ 7 


Neptune 


4.9818 x 10 6 


8 x 10~ 7 


Pluto 


3.66961 x 10 7 


3 x 10~ 6 



statistical uncertainties in the orbital elements, rescaled by a factor 10, to the time span of the 
fit performed by Pitjeva [13] . In view of future analyses, it is important to remark that eq. 
(PT7|) -eq, (fl8|) . assumed as real physical effects, show very distinctive patterns since there are 
no other known competing orbital effects on a and e, both Newtonian and Einsteinian. Indeed, 
general relativity does not predict any long-term rates of change for such orbital elements: the 
Schwarzschild field of a spherical static body affects just the pericenter and the mean anomaly 
of a test particle, while the Lense-Thirring effect due to the rotation of the central body causes 
secular precessions of the node and the pericenter. As far as Newtonian mechanics is concerned, 
the centrifugal oblateness of the primary does not impact a and e. Long-term rates of change 
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of the eccentricity may be induced by unmodeled/mismodeled distant distributions of matter 
like a massive ring or a pointlike object having certain specific orbital configurations, but their 
temporal signatures are known, and are quite different with respect to eq. (|18p so that it would 
be relatively simple to separate them. An isotropic mass loss suffered by the central body 
actually causes long-term variations of a and e, but, in the case of the Sun, they are smaller 
than the predicted cosmological values in Table [3l In any case, they have different patterns 
with respect to eq. (fT7]) -eq. ([IB]) . Moving to non-standard effects, in principle, a violation of 
the strong equivalence principle and spacetime variations of fundamental constants^ like, e.g., 
G can impact the eccentricity, but, again, not in the same way as predicted by eq. (|18p . 

In the following we will calculate the shifts per orbit Ay of various standard observables Y 

as 



2tt 



dY df dM dY dn 

^dM~df + ^~d^~dt 



where = 2'n/n^ is the orbital period along the osculating Keplerian ellipse, and M,I,u>, O 
are the mean anomaly, the inclination, the argument of pericenter and the longitude of the 
ascending node, respectively, of the orbit of the test particle. Y appearing in eq. (|2I [) is the 
analytical expression for the observable Y evaluated onto the unperturbed Keplerian orbit. 
dM/dt and dn/dt,K = a,e, I,u>, fi in eq. f)21 1) are the instantaneous rates of change of the 
osculating Keplerian orbital elements given by the right-hand-sides of the Gauss equations 
evaluated onto the unperturbed Keplerian osculating ellipse. Expressions for df /dM,dt/df can 
be found in standard textbooks. 

The shift per orbit of the radial distance 



2^ 



ail — e 

r = — — 

I + e cos / 

turns out to be 



(22) 



3iraH (l - e 2 ) (-2 + e 2 + 2y/l - e 2 ) o^„ rr ? 

Ar(™) = ) -±— 9 >- = JJ^E^ + o (e 4 ) . (23) 

n^e 1 4nb 

The radial velocity is left unaffected. 

The shift per orbit of the line-of-sight projection of the orbit of a binary system in the sky 

p = r sin /sin (u + /) , (24) 
which is practically determined from timing measurements of compact objects, is 

37raiT sini"sinu; (l - e 2 ) (-2 + e 2 + 2y/l - e 1 



^H 0S mI S mu f 1 _l\ +o{e4) (25) 



4re b V 2 



A time-dependent G(t) would change a as well, but not as predicted by eq. (|17[) . 
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where, in this case, I is the inclination of the orbital plane to the plane of the sky. Curiously, 
both eq. (f23|) and eq. ([25]) are independent of G because of eq. ([20]) . while they depend on 
M in such a way that they diverge in the limit M — > 0. Such a singularity might be explained 
by noticing that both eq. f)23[) and eq. ()25[) are the outcome of a perturbative calculation in 
which a Keplerian ellipse was assumed as unperturbed, reference orbit. Taking the limit M — > 
implies a breakdown of the validity of such an approximation since Ai H ^ cx v r oc 
for M — > 0, the Newtonian monopole = GM/r 2 becomes smaller than Ar . 
The shift per orbit of the radial velocity 



M; thus, 



n^a sin I [e cos ui + cos (oj + /)] 

Vp = • 

which is a typical spectroscopic observable in binaries studies, turns out to be 



2ttciHq sin / cos u 



2e 4 + 2 -1 + x/T 



+ e 2 -5 + 6Vl-e 2 



(26) 



biraeHo sin / cos cj 



i + ^) +0( ^). 



(27) 



Note that eq. (127ft is independent of M: in the limit M — > eq. (j27[) does not vanish. Instead, 
in the limit G — > eq. ([27]) correctly vanishes because of eq. ([201) . 



Theoretical motivations against the existence of a cosmolog- 
ical Hv r term in the equations of motion of test particles 



Actually, although appealing, the existence of eq. (fT6|) does not seem to be justified by a 
theoretical analysis of particle dynamics rooted in general relativity. 

The influence of the cosmic expansion on the gravitation fields surrounding individual ob- 
jects was the subject of several investigations since earlier times [III [T5[ fl6 [ [T71 [T8[ fl9 [ 1201 12T1 IB"1 
I22j ; for a recent review covering several aspects, see [3]. Connections with possible local varia- 
tions of cosmologically varying constants were elucidated by Shaw and Barrow in |23[ [2^1 [25| [26] . 
McVittie [H], working in the framework of general relativity, obtained a new spherically sym- 
metric metric describing a point mass embedded in an expanding spatially-flat universe from 
a suitable combination of the Schwarzschild and FLRW metrics. After a longstanding debate 
about its physical interpretation, it seems [271 [28] that the McVittie metric actually describes 
a point mass in an otherwise spatially-flat FLRW universe. However, Kopeikin [6], who quotes 
Carrera and Giulini [3], disagrees with such a conclusion. 

Indeed, let us consider the McVittie spacetime metric, written in "physical", i.e. not co- 
moving, coordinates {ct, r, 0,(ft} [3] [27] |2"2"] 



9oo 



1 - 2/i - Y 



(28) 



9oi 



h 



902 = 903 = 0, 



(29) 
(30) 



S 



where 



922 = ~r 2 , (32) 
g 33 = -r 2 sin 2 9, (33) 
912 = 913 = 923 = 0, (34) 



• GM , . 

A* = — 5~» 3 5 
err 



(30) 

c 



The equations of motion can be obtained from the Lagrangian 



C = - X -g^x\ (37) 



where the dot denotes derivation with respect to the proper time r, as 

d f dC\ dC 



0. (38) 



ar 



<0 



dr \dx^ J dxt 1 

As far as the radial equation of motion is concerned, the terms independent of c in dC/dr 



dC - 2 {—GM + H 2 r 3 ) . 2 Hri(l-3fi) m , 2 , 

= r^» 2 + ^ - >-t 2 - K -g + O c" 2 . 39 

dr H (i _ 2/i) 3/2 

The generalized momentum for r is 

• ^ r Hri 
Pr = ^ = —»-VT^- (40) 

In the weak-field and slow- motion approximation (1 3> /U,c — >■ oo,i — > 1), eq. (i39l) reduces 

to 

dC „ GM 



dr 



rc 



2 , tt2 



r 



2 



+ # 2 r - Hr, (41) 



while eq. (I40p becomes 

>r-Hr. (42) 

ar 

Thus, eq. yields 

f |^ j = r + g# 2 r + F 2 r - ifr. (43) 



2 As usual, we posed 6 = 7r/2. 



By equating eq. (03]) and eq. (JHJ) yields 



f - rcf 1 = — ( ^- - qH 2 r : (44) 

the left-hand-side is nothing but the radial acceleration in polar coordinates, while both the 
terms —Hr appearing in eq. (|43p and eq. (|4ip canceled each other in eq. (|44p . Arakida [29] 
recently studied the effects of the McVittie metric on the gravitational time delay. 

Another way to realize that standard general relativity does not predict the existence of 
eq. (|16p for a two-body system immersed in a FLRW expanding universe consists of looking 
at the cosmological impact on the binary system through the generalized Jacobi equation as a 
tidal effect in the local Fermi frame. [301 131] 132] . The generalized Jacobi equation is suitable 
for our purposes since it actually takes into account the relative velocity of the geodesies of 
the primary and the test particle freely moving in the background FLRW metric [20J . It turns 
out that an acceleration term linear in the velocity of the test particle with respect to the 
primary is, in general, present in the generalized Jacobi equation: it is proportional to [32| 
IZijkoV^ X k , i = 1,2,3, where X k and = X^ are the Fermi spatial coordinates of the local 
Fermi frame, and TZ^ upcT is the projection of the Riemann tensor of the background metric onto 
the Fermi frame. The problem is that the components of the Riemann tensor needed to yield 
eq. (|16p are vanishing in the FLRW metric [20] • It is well known that the standard Jacobi 
equation, with its TZjoko,j, k = 1, 2, 3 components, is able to yield the H 2 term of eq. (fT5]) |33j. 

Kopeikin [6], following a cosmological perturbation approach, did not find eq. (fl~6l) in the 
equations of motion of a test particle to the Newtonian level. He [6] found eq. (|15p . and a 
post-Newtonian term 

A (ff P N) =H ( v / c f v . (45) 

its magnitude is quite small in the Solar System. By employing the methods used in Section [3] 
it is possible to show that eq. ([4"5]) causes the following long-term variations of a and e 



a 



(HpN) 



2H a 3 n 2 (4 - 3y/l - e 2 j 2 H a 3 n 2 4H a 3 n 2 e 2 



c 



2 \/U? C 2 c 



+ + O (e 4 ) , (46) 



4if a 2 n?v / l — e 2 I — 1 + Vl — e 2 ) ou „2„2,, o- „2„2„3 
c z e c 2 2c 2 



Interestingly, eq. (j46|) -eq . (|47p do depend on both G and M in such a way that they correctly 
vanish in the limit G — )■ 0, M — )■ 0. About vanishing values of M, no singularity at all occurs in 
this case since ^4.( n )/^4(^p n ) oc Af~ 1//2 , so that the Newtonian monopole is always larger than 
the cosmological acceleration of eq. (|45p for M — > 0. Moreover, eq. (|46p is independent of the 
semimajor axis of the orbit of the test particle; to zero order in e, it amounts to the product 
of the Schwarzschild radius r s = 2GM/c 2 of the localized source times the Hubble parameter. 
For a Sun planet, eq. (|46p yields an increase rate of its semimajor axis as little as 21 jim. 
cty~ x + O (e 2 ). According to eq. ([47]) . the rate of change of the eccentricity depends on the 
size of the orbit as e/e ~ (r s /a)H. At first sight, the pericenter precession due to eq. ([T5]) . 
which amounts to [21] 

ui***) = _3go#o 2 v^? (48) 
2n b 
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gets singular in the limit G — > because of n h 1 in eq. (pl8|) ; actually, it is not so since [9] 



qH 2 = - = — p + 34 , (49) 



r9 . S A-kG ( p 
' -8—— (" + 3 S 

where p is the pressure of the total cosmic fluid, inclusive of any dark energy component. In 
principle, eq. (|48p has a singularity for M — > as well; actually, it does not make sense to 
consider such a limit since eq. (|48l) was obtained perturbatively [21] by assuming eq. (fl~5l) much 
smaller than the Newtonian monopole which, instead, would become smaller than eq. 

(|15p for M — > 0. The propagation of electromagnetic waves continuously exchanged between 
two atomic clocks geodesically moving in an expanding FLRW universe retains a cosmological 
imprint of order H large enough to allow for a possible detection in accurate range-rate Doppler 
experiments [BJ. According to Kopeikin [BJ, it may provide an explanation of the Pioneer 
anomaly [M] as well. 



5 General considerations on the viability of modified models of 
gravity from their observational consequences 

The discussions of the previous sections allow us to provide preliminary general criteria of 
viability of non-standard models of gravity [7] by inspecting the behaviour of their observable 
consequences. 

Let us suppose to have a given modified model of gravity 9Jt yielding a test particle extra- 
acceleration that can be viewed as a small correction to A^ in appropriate circumstances. 
Depending on the type of model, the parameter (s) entering A^\ collectively denoted as {a}, 
may or may not contain explicitly M and G. Let us, now, suppose to work out a certain observ- 
able effect X( m \ say an orbital extra-precession which slowly alters the otherwise unperturbed 
Keplerian ellipse. It is clear that, if we ideally switch off gravity, X^°^ must vanish in the limit 
G — > 0. Clearly, the same should occur in the limit M — > as well, unless certain conditions 
pertaining how is calculated must be satisfied. That poses certain basic constraints on 

the viability of the model and on the nature of its building blocks parameterized by {a}. It 
is important to note that it may happen that a well-behaved A^' yields a X®^ which, instead, 
is not. 

More precisely, let us consider the case of the Cosmological Constant A and the orbital 
effects it causes on a gravitationally bound binary system. The resulting particle acceleration 
is [35] 



= (50) 



which yields the orbital precession [35] [36] 



^ = ^nz. (5i) 

2n b 

Contrary to what might seem at first sight, eq. (|51jl is, actually, well-behaved with respect to 
G. Indeed, from [9] 

Ac 2 

« = 87rG (52) 

and eq. ([6j, it is possible to write 

Ac 2 = 3Q A H 2 . (53) 
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Since, according to eq. ([20]) . H 2 oc G in a spatially flat universe, we see from eq. (foTj) that 
o/ A ) oc y/G; thus, ci/ A ) — > in the limit G — > 0. The mass M of the central body enter eq. (|5ip 
as M -1 / 2 through the Keplerian mean motion nt>. As discussed in the previous Sections, it is 
meaningless to consider the limit M — > for the precession of eq. (|51[) since it was obtained 
perturbatively, while ^( N V-4 (A) < 1 for M 0. 

As another example, let us consider the Dvali-Gabadadze-Porrati (DGP) model [37]. It 
predicts an acceleration 

which is well-behaved with respect to the limits G — >■ 0, M — >■ 0; rg is a free length scale 
determined by observations of cosmological nature, while the =F sign depends on the cosmological 
expansions phases. As far as the observational consequences of eq. (|54p are concerned, it yields 
an orbital precession |38[ [39] 

. (DGP) 3C m i o N 

= T— + 0{c 2 ) (55) 

which is pathologically independent of both G and M. The same drawback is common to 
Galileon-based models as well [40 \ [4Tj yielding orbital precession [42] formally identical to eq. 

{55]). 

A further example is given by power-law extra-potentials [43[ 144] 

UW = ^,k> 1, (56) 

where R may or may not depend on GM along with other parameter (s) {a}. It can be shown 
perturbatively that the resulting orbital precessions are [45 \ [4*6] 

oc — ; (57) 
rib 

if ^ does not contain G, an unphysical singularity would appear in eq. (|57p for G — > 0. From 
this point of view, the Randall-Sundrum model [37] is well-behaved since & oc GM£ 2 in it, 
where £ is the anti-de Sitter (AdS) curvature scale. 

The Kehagias-Sfetsos (KS) |48] solution of the Hofava-Lifshitz (HL) modified gravity |49[I50] 
predict the existence of an additional radial acceleration 



4G 4 M 4 

where ipQ is a parameter of the HL model. The orbital precession induced by eq. (|58p is [51 



. , KS) 3G" 2 Af'/ 2 

it correctly goes to zero for G — > 0, M — > 0. The lack of singularity for M is due to the fact 
that the Newtonian monopole remains always larger than ^4( KS ) for vanishingly small values of 
M. 

Clearly, if a model is satisfactory from the point of view illustrated here, this does not 
necessarily imply it is really valid as a good description of (a part of) the physical reality: it 
has to pass other independent tests. 
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6 Conclusions 



We investigated the possibility that the cosmological expansion may impact the orbital motion 
of a localized two-body system at first order in the Hubble parameter H. 

Reasoning classically, intuitive guesses lead us to postulate, at Newtonian level, the existence 
of an additional radial acceleration Ai HN ^ of order H proportional to the radial velocity v 
of the proper motion of the test particle with respect to the primary. By considering ^4^ N ^ as a 
small correction to the Newtonian monopole A^ N \ we perturbatively worked out its long-term 
effects on the osculating Keplerian orbital elements of the test particle and on its distance from 
the primary, and, in the case of a binary in the plane of sky, on the projection of its orbit 
onto the line-of-sight and on its radial velocity. It turned out that both the semimajor axis a 
and the eccentricity e of the test particle would secularly increase. The analytical expressions 
of their rates of change are independent of the mass M of the primary, so that they do not 
vanish in the limit M — > 0, contrary to that one might reasonably expect. Such a feature 
constitutes a general requisite to be satisfied in the search of potentially viable alternative, i.e. 
non-cosmological, physical mechanisms able to provide Ai HN \ 

Nonetheless, the consequences of Ai HN ^ would be interesting from a phenomenological point 
of view since their magnitude is close to the current level of accuracy in determining the plane- 
tary orbits in our Solar System. Preliminary analyses of recently published data show a tension 
among the predicted effects for Ai HN ^ and the admissible upper bounds on the secular rates of 
change of a and e for Mercury and Mars. 

Then, we looked for theoretical justifications of the existence of Ai H ^ in three different 
frameworks within general relativity. None of them provided it at Newtonian level, contrary 
to the known Hooke-like term of order H 2 . Instead, as recently pointed out in literature, a 
velocity-dependent acceleration ^4^ pN ) of order H and directed along the velocity v of the 
test particle exists at post-Newtonian level. We perturbatively worked out its orbital effects 
by finding secular rates of change of a and e proportional to the Schwarschild radius r s of the 
primary. They are well-behaved since they correctly vanish in the limit M — > 0. For a planet 
of our Sun such effects are negligibly small: suffice it to say that the semimajor axis increases 
at a rate of just 20 /j,m per century. We discussed the limits of validity of the orbital precession 
due to the H 2 term. We noticed that the formal singularity occurring in it for M — > actually 
has no physical meaning since such a limit is beyond the regime of validity of the perturbative 
calculation yielding the precession itself. 

Finally, we extended some points emerged in dealing with the previous cosmological issues 
by discussing certain basic criteria of viability that modified models of gravity should 
generally meet in view of their predicted observable effects like, e.g., orbital precessions. 
Given that such modified gravities must reduce to small perturbations of standard Newtonian 
gravity in appropriate circumstances, their precessions, which slowly alter an otherwise fixed 
Keplerian ellipse, must necessarily vanish in the limit of no gravity, i.e. for G — > 0. The same 
should occur also for M — > 0, unless such a limit violates the validity of the pertubative 
regime in which the precessions are calculated. If a given model provided a precession which 
is, say, independent of G and/or M, caution should be used in further considering that model. 
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